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ABSTRACT 



Context, While the sources of X-ray and radio emission in the different states of low-mass X-ray binaries are relatively well under- 
stood, the origin of the near-infrared (NIR) and optical emission is more often debated. It is likely that the NIR/optical flux originates 
from an amalgam of different emission regions, because it occurs at the intersecting wavelengths of multiple processes. 
Aims. We aim to identify the NIR/optical emission region(s) of one such low-mass X-ray binary and black hole candidate, 
XTE J1650-500, via photometric, timing, and spectral analyses. 

Methods. We present unique NIR/optical images and spectra, obtained with the ESO-New Technology Telescope, during the peak 
of the 2001 outburst of XTE J1650-500. 

Results. The data suggest that the NIR/optical flux is due to a combination of emission mechanisms including a significant contribu- 
tion from X-ray reprocessing and, at early times in the hard state, a relativistic jet that is NIR/radio dim compared to similar sources. 
Conclusions. The jet of XTE J1650-500 is relatively weak compared to that of other black hole low-mass X-ray binaries, possibly 
because we observe as it is being "turned off" or quenched at the state transition. While there are several outliers to the radio-X-ray 
correlation of the hard state of low-mass X-ray binaries, XTEJ1650-500 is the first example of an outlier to the NIR/optical-X-ray 
correlation. 
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1. Introduction 

Transient low-mass X-ray binaries (LMXBs) are, for most of 
the time, in a state of quiescence with faint or non-detected 
X-ray emission and near-infrared (NIR)/optical emission dom- 
inated by the main-sequence companion star (possibly with sig- 
nificant contribution from the cold accretion disk). They are of- 
ten only discovered when — powered by an increased level of 
accretion onto the central compact object (black hole or neutron 
star) — there is a dramatic increase of the X-ray, NIR/optical 
and radio flux. During these outbursts the systems are observed 
to go through several high-energy (X-ray) spectral states before 
returning to a quiescent state, usually on timescales of weeks, 
months, or even longer. 

LMXBs with a black hole compact object are observed (in 
X-rays) initially in a generally low-intensity, power-law dom- 
inated, hard state before transitioning to a usually higher in- 
tensity, thermal- dominant, soft state that decreases in flux and 
evolves via a late hard states back into a quiescent state. Some 
outbursts also display a steep power law state with a steeper 
spectral slope than the regular hard state during the transition 
between hard and soft states. In general, the X-ray hardness ver- 
sus intensity diagra ms (HIDs) are observed to follow canonical 
trajectories (e.g., iHoman et al.l l200lt iHoman & Bellonil [2005; 



* Based on observations made with the European Southern 
Observatory telescopes obtained from the ESO/ST-ECF Science 
Archive Facility. 



Bellonil I2010h and are often (though not always) indicative 
of the black hole nature of the compact object. Additionally, 
the hard states are associated with aperiodic variability of the 
light curve not present in the soft state and, in hard and steep 
power law states, quasi-periodic oscillations (QPOs) are de- 
tected (for a fuller description of the various possible states see 
iMcClintock & Rermllardl2006l) . 

The origin of the NIR, optical and ultraviolet (UV) emis- 
sion from a black hole LMXB is much less understood than the 
origins of the X-ray or radio emission because the optical wave- 
lengths are at the intersection of a number of different emission 
mech anisms (for reviews of optical pr o perties of LMXBs se e 
e.g., Ivan Paradiis & McClintockl Il995t (Charles & Coel l2006h . 
The companion star is generally relatively dim and may not 
contribute significant flux during outburst when other mecha- 
nisms are active. Intrinsic, thermal emission from the hot, outer 
accre tion disk may contribute at UV and optical wavelengths 
(e.g., IShakura & Sunvaevlll973l iFrank et alJ [2002h though the 
reprocessing of X-rays in the same region of the accretion 
disk is thought to be a significant sourc e of flux at wave- 
lengths from UV through to N IR (e.g., ICunninghaml 1 1976b 
Ivan Paradiis & McChntockl [l994). Recently, evidence has been 
mounting that the relativistic jet, usually detected in radio, also 
produces a significant contribution to the NIR flux, at least 
in the hard state (e.g.. iJain et alJl200lb ICorbel & Fended l2002t 
IChatvet all 120031: iRussell et al.l I2OO6I) . and it is possible that 
the power law X-ray emission in the hard state extends to, and 
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Table 1. Nights of observations 



Night 


MJD 


Filters 


Size T (') 


ESOID 


2001-09-08 


52161 


J,H,K S 


5x5 


67.B-0486 






GBF, GRF* 


— 


67.B-0486 


2001-09-14 


52167 


V,R,l 


9x9 


67.B-0315 


2001-09-15 


52168 


RILD#l* 




67.B-0486 


2001-09-24 


52177 


J,H,K S 


4x4 


67.D-0200 


2001-10-04 


52187 


J,H,K S 


5x5 


68.D-0316 


2001-10-06 


52189 


K s 


4x4 


68.D-0316 


2001-10-09 


52192 


V,R,I 


9x9 


68.A-0440 


2001-10-10 


52193 


RILD#l* 




68.D-0144 


2001-10-13 


52196 


V 


1 x 1 


68.D-0144 


2003-05-13 


52772 


I 


3x 10 


71.D-0337 



1 Final image size after cropping; * Spectra. 



contributes at, optical wavelengths. This emission is ascribed to 
the inner regions of the system — either a corona, advection- 
dominated accretion flow, or the base of a compact jet — which 
we will refer to as the corona for simplicity. 

The black hole X-ray transient XTE J 1650-500 was ini- 
tially detected by the Ro ssi X-ray T iming Explorer (RXTE) 
on September 5 2001 dRemillardl 1200 ll) . In the follow- 
d ays, its variable nature was confirmed 
Markwardt et afl 120011: IWiinands et al., 



Remillard, 



oscillations (QPOs) indicative 



200 lh 
of black 



mg 
l200ll 

quasi-periodic 

J evnivtsev & Sunya evj 12001 
OOlt lAueusteiin et all [200 it 
iBuxton & LondishU200it) and ATCA radio (iGroot et all 12001) 
counterparts were also reported at the time. The canonical black 
hole state behaviour of the source over th e ~ 100 day Ion 
outburst i s well documented in both X-ray dRossi et all 1200' 
and radio (ICorbel et all 120041) . while evidence of the black hole 
nature of the compact o bject is offered by the confirmation of 
QPOs dHoman et alll2003b. the detection of a rel ativistic Fe Ka 
line (e.g.. iMiller et al.ll2002t iMiniutti et al.ll2004l) and a derived 
mass of ^ 5.1Mm dSlanv & Stuchlfldl2008h . 

The optical observations of the source during outburst (cited 
above) wer e limited to a few epochs , but o n June 10 2002, in 
quiescence. lSanchez-Fernandez et al] f2002) obtained 15 optical 
spectra with FORS2 on the 8.2m VLT at ESO, Paranal. Using 
these spectra and relative optical photometry (obtained with the 
6.5m Clay telescope at Las Campanas Observatory, Ata cama 
during May, June, and August 2002), lOrosz et all d2004l) de- 
rive an orbital period of 0.32 days (7.69 hours). On August 
2 2002 iGarcia & Wilkes! d2002l) estimated the quiescent opti- 
cal magnitude. In this paper we present the only significant 
NIR/optical observations during outburst, obtained by the ESO 
NTT in September and October 2001 (TableQ]), which comprises 
all available unpublished, archived ESO data of the source. In 
section [2] we introduce the observations and reduction methods, 
while in section[3]we present the results of our photometric, tim- 
ing, and spectral analyses of the data. We discuss the interpreta- 
tion of our findings within the context of black hole LMXBs in 
section|4]and summarise in section|5] 

2. Observations and reduction 

2.1. Photometry 

Optical (V, R, I) and NIR (J, H, K s ) dat a were obtained wi th the 
ESO Multi-Mode Instrument (EMMI: [Dekker et alj[l986h and 
the Son of ISAAC (Soft) infrared spectrograph and imaging 




Fig. 1. NTT 60" x 40" finding charts (upper: 970s V image; 
lower: 910s Ks image) with the 2" optical positional uncertainty 
dCastro-Tirado et ail 1200 lh marked by a circle. 



camera on the 3.58m ESO - New Technology Telescope (NTT), 
respectively, on a number of nights in September and October 
2001, and in May 2003 (Table[TJ. Data were reduced using the 
IRAF package, in which crosstalk correction, bias-subtraction, 
flatfielding, sky subtraction, bad pixel correction and frame ad- 
dition were carried out as necessary. NIR data obtained on the 
nights of September 8 and October 4 did not use a dithered point- 
ing pattern as is usual for NIR images, therefore it was not pos- 
sible to produce a reliable sky for subtraction; instead the sky 
produced on September 24 was used. This should be of min- 
imum effect because we used relative photometry to estimate 
magnitudes. 

The images were as trometrically calibra ted against 2MASS 
dSkrutskie et all [2006b or USNO-B1.0 dMonet et all [2003h 
within the GAIA package. The position of the source was de- 
rived via the point spread function (PSF) photometry of the deep 
Ks band image on October 6 (MJD 52188; seeing as 1.2") as 
16:50:00.95 -49:57:44.34 with a positional erroiQ dominated by 
the 0.1" 2MASS systematic uncertainty (Figure[T]i. 

Relative PSF photometry was carried out on the final im- 
ages using the DAOPHOT package (IStetsonl 119871) within IRAF. 
However, the optical data obtained on September 14 and October 
9 were so poorly binned that producing an acceptable PSF model 
was impossible (full width at half maximum, fwhm < 2 pix- 
els); hence these images use relative aperture photometry. The 
NIR magnitudes (Table[2j Figure[2]i were calibrated against the 



All uncertainties in this paper are given with a confidence of lcr. 
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Table 2. Optical and NIR exposures and magnitudes 



A 


*5 


* 


H 





J 


■ 


I 


• 


R 


o 


V 



Hard 



o 



Inter. 



Steep PL 



XI 
Soft 



2160 



2170 



2180 



2190 



MJD 



Fig. 2. V, R, I, J, H, and Ks band light curves (as presented in 
Table |2]i with the pe riods of the hard, in termediate, steep power 
law and soft states (ICorbel et all 2004) marked. \A\ marks the 
epochs of the spectral observations. 



2MASS catalogue using « 200 objects per image, after outliers 
were removed. The calibration, at least for the Ks band, was 
confir med on the October 6 Ks image and four iPersson et al.l 
(1998) photometric standards, observed on the same night (the 
only night on which reliable NIR standards were available). 

Optical magnitudes (Table|2j Figure|2]i were calculated rela- 
tive to «10 relatively isolated field stars, which we re calibrated 
on the September 14 images against Landolt d 19921) photometric 
standards, observed on the same night (the only night on which 
optical standards were obtained). The errors on the optical mag- 
nitudes are dominated by this absolute photometric calibration 
but unfortunately, due to the saturation of USNO-B 1 .0 objects 
in the field, we are unable to calibrate against that, or other, cat- 
alogues. The upper limit of the late / band image is approxi- 
mated from the dimmest observable object in the region of inter- 
est. Because we were only able to obtain aperture photometry on 
September 14 and October 9, the measured magnitude is likely 
to suffer contamina tion from the unresolved, nearby sources (see 
lOrosz et all d2004l) for a better resolved image from the 6.5m 
Clay telescope). To correct for this, we estimated the magnitude 
of the contaminating sources to be / = 19.2 + 0.5 from the mag- 
nitude at the position in the late / band image (May 13 2003) 
when the source had faded; V = 18.4 + 0.4 from the residual of 
the PSF photometry on October 1 3 200 1 ; R = 1 8 .7 + 0.4 from an 
interpolation of the / and V contaminating magnitudes. The dif- 
ference in magnitude ranges from > 1 magnitude for the V band 
to >3 magnitudes for the / band. In each case, given the rela- 
tively large absolute calibration error, the corrected magnitude 
(Tabled is consistent, within errors, with the magnitude uncor- 
rected for contamination. 

The data obtained on October 6 and October 13 consisted of 
multiple high-cadence, "fast" photometry images, in Ks and V 
bands, respectively, to investigate possible short-term variability 
(Tabled up to ~ 2-3 hours, with an approximate sampling of 
an image every 90s. PSF photometry was carried out on each of 
these individual images, again using DAOPHOT. Source magni- 
tudes were calculated relative to a number of field stars (10 in 
NIR and 8 in optical) and normalized so that the weighted aver- 
age is equal to zero. In each, a few magnitudes were eliminated 



MJD-50000 


Filter 


N x Exp (s) 


Magnitude 


2167.00354 


V 


1 x30 


17.1 +0.2 


2191.99546 


V 


1 x 120 


17.3+0.3 


2195.99182 


V 


97 x 10 


17.4 + 0.4 


2166.99704 


R 


1 x 150 


17.0 + 0.4 


2167.00545 


R 


1 x30 


16.7 + 0.4 


2191.99848 


R 


1 x 120 


17.0 + 0.4 


2192.00531 


R 


1 x40 


17.3+0.4 


2167.00002 


1 


1 x60 


15.86 ±0.13 


2167.00699 


I 


1 x30 


15.69 ±0.13 


2192.00104 


I 


1 x90 


16.3 ±0.3 


2192.00696 


I 


1 x40 


16.0 ±0.2 


2772.17920 


I 


4x300 


>19.0 


2161.02124 


J 


5 x20 


14.42 ±0.11 


2176.98651 


.1 


9x 10 


14.65 ± 0.09 


2186.98230 


J 


5 x20 


14.79 ±0.10 


2161.02907 


H 


5 x20 


13.79 ±0.11 


2176.99397 


H 


9x 10 


14.16 ±0.12 


2186.99013 


H 


5 x20 


14.12 ±0.16 


2161.03691 


K s 


5 x 20 


13.29 ±0.13 


2177.00145 


Ks 


9x 10 


13.82 ±0.12 


2186.99797 


Ks 


5 x20 


13.99 ±0.16 


2188.97400 


K s 


130x7 


14.18 ±0.13 


2160.98738 


GBF 


2x 180 




2160.99906 


GRF 


2x360 




2167.99490 


RILD#\ 


2x600 




2193.01287 


RILD#\ 


2x900 





due to unsatisfactory PSF subtraction residuals which became 
obvious during visual inspection of the images. 



2.2. Spectral energy distributions 

The observed magnitudes (Tabled or, if there were multiple 
observations per band per night, the weighted average magni- 
tudes, were first converted into flux densities, F v , at frequency 
v. They were then converted into flux per filter, Ffu ter in units 
of photons crrT 2 s . This conversion was made via Ffn ter = 
1509.18896F V (AA/A) where A and AA are the effective wave- 
length and width of the filter in question. XSPEC compatible files 
for spectral energy distribution (SED) fitting were then produced 
from the flux per filter value using the FTOOL, f lx2xsp. In fig- 
ure [3] we c orrect for the Galactic extinction in the direction of 
the sourc e dSchlegel et all 1 19981) of E B - V = 0.923 (A K ~ 0.3, 
A v ~ 2.9; Car delli etalj fl989fl This value of extinction is con- 



sistent with that implied (Guver & Ozel, 2009) from the X-ray 
absorption of this source (Montanariet al., 2009), even though 
the source is likely closer dHoman et all 120061) and should not 
suffer the full effect of Galactic dust. 



2.3. Spectroscopy 

Optical and NIR spectra were obtained with EMMI and SofI 
on nights in September and October 2001 (Tables[TJ[2]i. EMMI 
obtained red (3, 850-10,000 A) low-dispersion spectra using 
grism#l (RILD1) on two nights. SofI obtained usable blue 
(GBF; 9,500-16,400A) and red (GRF; 15, 3 00-25, 200 A) low- 
resolution spectra on a single night, while spectra from another 
two nights (September 18 and October 4 2001) were not use- 
ful due to a very high level of spatial dispersion, as well as a 



2 This extinction should be treated with caution as estimates so close 
to the Galactic plane (< 5 deg) are uncertain. 
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Fig. 3. Flux versus frequency plot at five different epochs, uncor- 
rected (upper) and corrected (lower) for the Galactic extinction 
in the direction of the source of Eg-v - 0.923. 



lack of spectra for wavelength calibration on those nights. The 
data were reduced using the IRAF package in which crosstalk 
correction, flatfielding and dark- or bias subtraction were carried 
out as necessary. To correct for NIR sky, the dithered NIR expo- 
sures were subtracted from each other and the resulting spectra 
summed later. 

Spectra were extracted and reduced within the IRAF pack- 
age, noao . twodspec. Wavelength calibrations were performed 
against helium + argon (optical) or xenon (NIR) lamps whose 
spectra were extracted using the same parameters as for the rel- 
evant source. Atmospheric telluric features significantly affect 
the spectra and were corrected for by dividing the source spec- 
trum by that of a telluric standard at a similar airmass, using the 
telluric tool within IRAF. This procedure often causes arte- 
facts in the corrected spectra and, in the case of our optical spec- 
tra, these artefacts are significant, which is why we did not apply 
this procedure. Neither the optical nor the NIR spectra were flux 
calibrated, but they were normalized. 

The NIR spectra (from the night of September 8 2001) were 
combined with the optical spectrum from the night of September 
15 2001 and the final spectrum, from 3, 850-25,200 A, is plotted 
in Figure [4] Owing to discrepancies in the optical wavelength 
calibrations between the two nights, we are unable to sum the 
two optical spectra and instead plot the spectra with the lower 
noise. The (wavelength-dependent) resolution of the final spec- 
trum is 7-10 A, with a wavelength calibration error of < 20 A 



(optical) or < 40 A (NIR). The apparent features just redward of 
4000 A and blueward of 10 4 A are at the edge of the optical grism 
range where the response is not reliable. The high levels of noise 
at ~ 14, 000 A and ~ 19, 000 A are artefacts of the correction for 
telluric features, while the general trend towards higher noise at 
the red end is a real effect of the observations. An absorption 
feature with a fwhm of 60A at 7,625 A is detected in the optical 
spectrum of the night of September 15, but this is an atmospheric 
feature, also detected with similar properties in the spectroscopic 
standard on that night. This feature is not detected in the October 
10 spectrum, though this spectrum has a root mean square noise 
about four times greater than the September 15 spectrum. 

3. Results 

3.1. High-cadence photometry 

Though the high-cadence light curves (Figures [5] and [§}, partic- 
ularly the Ks band, seem to display some level of variability by 
eye, this is not supported by a careful analysis of the data. The 
standard deviation of the V data is 0.012 magnitudes, compared 
to an average error on the data points of 0.01 1 ± 0.002. We fit- 
ted a constant value to the light curve which, for the V band, 
returns xl = 143.1/(97 - 1) = 1.49, statistically inconsistent 
with being an acceptable fit at >99% confidence. However, if 
we apply the same tests to a number of nearby sources of sim- 
ilar magnitude (and hence similar signal-to-noise ratio) we find 
similar standard deviations (0.01 1 - 0.014) and^j; (1.31 - 2.05), 
for the same number of degrees of freedom. This implies that 
any variability is not intrinsic to the source itself but is likely 
noise. Likewise, the Ks band data, which have a standard devi- 
ation of 0.059 magnitudes compared to an average error on the 
data points of 0.017 + 0.007, returns^ = 1222/(128- 1) = 9.6, 
inconsistent with a non-variable source. However,several nearby 
sources of similar magnitude also return similar standard devi- 
ations (0.056 - 0.086) and^ (9.1 - 14.9), again implying that 
any variability is not intrinsic to the source itself and that, for the 
Ks band, the errors on the individual photometric points are un- 
derestimated, likely due to e.g. the difficulties of sky subtraction 
and PSF modelling in such a crowded field (FigureQ]). 

Furthermore, we used the IRAF task, pdm — an imple- 
mentation o f the phase dispersion minimization method of 
IStellingwerj (|J978) — to test for periodicity in the variability 
of the light curves (see periodograms in lower panels of figures 
|5]and|6]). We find no reliable periodicity in either of the observed 
bands as the tests returned Stellingwerf statistics, 0^1 for all 
periods less than the length of observations and broad minima 
of > 0.7 for longer test periods. While the statistic in each 
band was generally lower than that for the nearby comparison 
objects (introduced above, dashed lines in periodograms), there 
was no strong local minimum that would sugest a periodic vari- 
ability. Neither is there any sign of periodicity when the statis- 
tic of the source is normalized relative to that of the comparison 
objects. This null result was also confirmed using the Lomb- 
Scargle method for unevenly sampled data dPress & RvbickiL 
1989), which shows no significant peak in power over the pe- 
riod ranges in question. 

We also tested both light curves to see if either were consis- 
tent w ith a sinusoidal periodicity similar to that of lOrosz et al.l 
(2004), who measured a quiescent period of 0.3205 days (7.69 
hours, compared to the ~2-3 hour duration of our observations) 
and an R band peak-to-peak amplitude of 0.2 magnitudes. This 
fit to the Ks light curve returns^ = 7.2 (125 degrees of free- 
dom) while similar fits to the comparison objects return xl of 
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5,000 
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20,000 



25,000 



Fig. 4. Final normalized spectrum from 3,850 to 25,200 A. The apparent features just redward of 4,000A and blueward of 10, 000 A 
are at the edge of the optical grism range where the response is not reliable. The high levels of noise at ~ 14, 000 A and ~ 19, 000 A 
are artefacts of the telluric correction while the absorption feature at 7,625 A is atmospheric. 



between 10.3 and 14.2. A similar fit to the V band light curve 
of the source returns xl - 1-5 (93 degrees of freedom) while 
fits to the comparison objects return^ of between 1.5 and 3.4. 
While it would seem that both bands could be consistent with 
this period, it is necessary in both cases that the epochs of obser- 
vations fall serendipitously at the minima of the periodic func- 
tion, which seems unlikely. Indeed, to observe such a period 
in outburst would be physically unrealistic since the periodicity 
was detected through orbital eclipsing; during outburst, where 
the optical emissio n is 5 magnitudes brigh ter than quiescence 
(V ~ 24, R ~ 22: iGarcia & Wilkes! 12002). this would not be 
detected. However, so-called "superhump" modulations, indica- 
tive of the orbital period, ha ve been observed in the hard state of 
several X-ray binaries (e.g JZurita et al.l f2008). 



3.2. Spectral energy distributions 



Though the light curves (Figure|2]i and fiux-versus -frequency 
plots (FigureO display apparent hints of NIR colour evolution 
as well as a possible V band excess, the spectral indices (a, 
where F v oc V) of the two sets of V,R,I data and the three 
sets of J, H, Ks data are all consistent within the (not particu- 
larly well constrained) errors. Furthermore, a simultaneous fit of 
all five (non-simultaneous) SED epochs does not support signif- 
icant spectral evolution, as the data are well fit (xl = 0.72 for 
8 degrees of freedom) by an absorbed power law of spectral in- 
dex, a = 0.8 + 0.2, with the extin ction set to Eb-v = .923, 
the Galactic value in that direction (Schleg el et all Il998h . The 
only evolution, at least given the quality of the data, is that of the 
flux offset, which decreases in time, as expected from the light 
curves. It is important to note here that while spectral evolution 
is not statistically supported, neither is it ruled out as a possibil- 
ity. 



Alternatively, the combined SEDs could also be acceptably 
fit by simple black body radiation at any temperature > 3100K 
(> 2.7x 10~ 4 keV) with extinctions, E B _ V < 1 .6, or at 7000 K (6x 
10 _4 keV) for the Galactic extinction (xl = 1.67). Additionally, 
the SEDs are also consistent with black body radiation from an 
accretion disk of temperature > 5500K (> 4.7 x 10~ 4 keV) with 
extinctions, E B - V < 0.8, or at a; 6 x 10 6 K (0.5 keV) for the 
Galactic extinction (xl = 1.86). It is hence clear that the SEDs 
are not particularly constraining and results should be treated 
with caution, especially because of the time-averaged nature of 
the above fits. 



■a 




0.5 1 
Time (hours) 



'5 o 




2 4 6 

Period (hours) 

Fig. 5. Top: High-cadence V band light curve of the source along 
with three comparison objects of a similar magnitude (normal- 
ized to a relative magnitude of 0. 1) and, bottom, the periodogram 
(© statistic versus period) for the source (solid line) and the 
nearby comparison objects (dashed lines). 



3.3. Spectroscopy 

The final spectrum (Figure |4|i displays a single possible fea- 
ture that we might attribute to the source: a tenuous line (with 
a peak at ~ 6 times the noise and afwhm of 10 A) is detected 
at 16,700 A. The same feature is visible in the spectrum before 
telluric correction, therefore is not an artefact of that process. 
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Fig. 6. Top: High-cadence K$ band light curve of the source 
along with three comparison objects of a similar magnitude (nor- 
malised to a relative magnitude of 0.6) and, bottom, the peri- 
odogram (0 statistic versus period) for the source (solid line) 
and the nearby comparison objects (dashed lines). 

This line could correspond to an HI line at 16,8 10A (Brackett 
series, rj), though given the difference in wavelength and the 
very narrow fwhm compared to the spectral resolution (7-10 A), 
it is likely that this line is not real. If this line were real we 
would also expect to detect other lower energy Brackett series 
lines that would have higher abundances, but we do not de- 
tect any , again suggesting that th e line may be a statistical fluc- 
tuation. lAugusteiin et alJ (120011) . in their optical spectrum ob- 
tained on September 8 2001 (on the night of our NIR spectra 
and seven nights prior to our low-noise optical spectrum) identi- 
fied Balmer series emission lines, but we observe no such lines 
with a signal-to-noise ratio > 4. Neither do we detect the broad, 
double-peaked Ha emission line from the accretion disk or the 
weak absorption lines from th e companion star identified by 
Sanc hez-Fernandez et alJ d2002l) in their VLT spectra from June 
10 2002 (i.e., in quiescence). 

4. Discussion 

The states of XTEJ1650-500 are defined, from X-ray and 
radio observations, as being hard up to MJD -52,161; 
hard/intermediate up to MJD ~52,172; steep power law up to 



MJD -52,187; ther mal-dominant/soft u p to MJD -52,232; and 
hard after that date (ICorbel et all 12004). Our observations were 
hence obtained in the very late, initial hard/intermediate, steep 
power law and in the first few days of the soft state (Figure|2]i. 
Specifically, our multi-band NIR observat ions were ob tained 
within 5 hours of observations 2, 3, and 4 of lCorbel et all in the 
hard and steep power law states. In figure|7]the dereddened NIR 
and radio data are plotted along with the X-ray flux from that 
paper (converted to Jansky at 4.7 keV). The common spectral in- 
dex of the NIR observations (section USl i is also plotted along 
with the fitted spectral indices of the radio (ICorbel et al.l 120041 
who used a therma l free-free absorbed power law) and the X-ray 
(Rossi et al., 2004) at those epochs. One can see that epoch 1, in 
the hard state, is characterised by a high level of radio emission 
and a relatively shallow-power-law dominated X-ray flux, while 
epochs 2 and 3 display declining radio flux and steep power- 
law-dominated X-ray emission. One of the most obvious aspects 
of this plot is that while the radio flux drops by over an order 
of magnitude, the X-ray flux remain s relatively st able, dropping 
only fractionally. As pointed out bv lCorbel et all this is consis- 
tent with the radio originating from a self-absorbed compact jet 
while also indicating that the jet contributes little, if any, flux to 
the X-ray emission. Furthermore, the decrease to background of 
radio flux is consistent with the assumed behaviour of the jet in a 
black hole system that is quenched duri ng the soft, thermal state 
( ICorbel et alll2004tlFender et alll2004l) . 

4. 1 . The source of NIR/optical emission 

It is expected that the NIR/optical emission of a black hole 
LMXB originates from a combination of the companion star, the 
relativistic jet, intrinsic thermal or reprocessed X-rays from the 
accretion disk, and a possible contribution from the corona (see 
section [TJ. The estimate d quiescent optical m agnitudes of this 
system (V ~24,R~ 22 lGarcia & Wilkesll2002t / > 19.0 herein) 
are significantly (-100-500 times) dimmer than those observed 
in outburst and hence imply that the companion star has no mea- 
surable effect on the emission during outburst. If the observed 
NIR flux were due exclusively to the jet, we would expect an 
almost flat spectral index (a ss 0), much shallower than that ob- 
served (a = 0.8 + 0.2). We would also expect an approximately 
linear relationship between the NIR and radio flux (which is as- 
sumed to be, and is consistent with being, due only to jet emis- 
sion). However, it is unclear how this relationship would evolve 
over a state transition, as is being observed here, since flux at one 
frequency may decrease earlier or faster than that at another fre- 
quency. We find that the NIR flux falls off much more gradually 
with Fnir 00 ^RaJio 07 ( at ^ otn ra di° frequencies), and maintains 
a high level of emission even as the jet/radio emission goes to 
zero in the soft state. Additionally, over the same time, the X- 
ray flux falls off even less steeply than the NIR (Fnir k F ^ 5±1 °), 
therefore it seems plausible that at least some of the NIR emis- 
sion originates from the jet while the remaining emission must 
be due to the accretion disk, reprocessing and/or coronal emis- 
sion. 

Unfortunately the X-ray flux cannot be extrapolated to imply 
an accretion disk contribution to the NIR, as the accretion disk is 
not generally observed to make a significant contribution to the 
X-ray in the hard or steep power law states. We would expect to 
observe spectral emission lines if the NIR emission originated 
from an accretion disk, but we observe none in our spectra dur- 
ing the hard/intermediate states; however, this is not to say that 
the accretion disk is not making a, possibly significant, contribu- 
tion at these times. Neither can we extrapolate the X-ray flux to 
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Fig. 7. Dereddened/unabsorbed (except for radio) fiux-versus-frequency plot at three epochs (1 = MJD 52161, black; 2 = MJD 
52177, red; 3 = MJD 52187, blue). Spectral indices, a, are from previously published sources, except for the NIR. Power laws are 
purely phenomenological in order to show the relative behavior at different frequencies. 



imply the corona's contribution because, while it is clear that the 
power law contribution in the hard and steep power law states 
must break at low energies, we cannot constrain the form of this 
break. Since the significant spectral change of the X-ray does 
not seem to affect the NIR spectra, which remain constant, one 
might be tempted to imply that the corona does not contribute to 
the NIR, but because of the poorly constrained nature of the SED 
and the lack of knowledge of the low-energy spectral shape of 
the corona, this is not necessarily the case. As mentioned earlier 
(section [1), black hole LMXBs display significant X-ray vari- 
ability in the hard state at frequencies from ~ 10~ 3 to 10 4 Hz, 
which is absent from the soft state, and this hard state variabil- 
ity h as also been observed in the optical fro m ~ 10~ 4 to 10 Hz 
(e.g.. iCasella etal]l201(* IChatv et id] 1201 lh . The high level of 
variability in both optical and X-ray is often attributed to the 
non-thermal emission of the jet or corona. We find no such vari- 
ability (section [Til l in our high-cadence optical and NIR obser- 
vations (spanning ~ 10~ 4 to 10 _1 Hz), which might be indicative 
of emission associated with these regions, consistent with the 
source's soft state at those epochs. 

If the NIR emission was primarily due to X-ray reprocessing, 
the observed NIR flux s hould be proportional to the ob served 
X-ray flux as oc F~° 5 (Ivan Paradiis &McChntockll 19941 or for 
more detailed calculations. ICoriat et al. 2009) but instead it fol- 
lows the much steeper relationship of Fmr k F^ ±w . However, 
it must be noted here that the level of reprocessing (and the 
NIR - X-ray relationship) is very sensitive to the X-ray spec- 
tral shape, and hence the X-ray range sampled by a particular 
instrument. Over a state transition, as is being observed here, 
the X-ray spectrum is undergoing dramatic changes of spectral 
shape; consequently, while the X-ray flux over the measured 
range remains constant, the flux over the range contributing to 



the reprocessing may be significantly variable. We do find that 
the absolute optical magnitude of the source (My = 2.2) assum- 
ing a distance of 2.6 + 0.7 kpc (Homan et alll2006h . agrees well 
with the observed correlation with £ = (L x /L Edd y /2 P 2/3 =0.4 
(Ivan Paradiis & McClintocl Il994t iDeutsch et all |2000|) . when 
we use the observed period, P = 7.69 days (lOrosz et al] 
and the derived mass of ^ 5.1M Q dSlanv & Stuchllk, 



2004) 



2008) 



This correlation was derived for reprocessing, implying that the 
optical flux may originate from this process. Furthermore, the 
lack of observed variability in the long term optical — in con- 
trast to the NIR — light curves (Figure|2]i is consistent with 
the ^optical x F x 5 relationship expected from reprocessing. 
The possible, though not statistically significant, V band excess 
(Figure[3]l may be due to reprocessing peaking in that, or a bluer, 
band. 



4.2. The weak jet of XTE J1 650-500 



It has been noted dCorbel et al.L l2004t IFender et all [2010) that 
XTE J 1650-500 is X-ray loud, relative to th e radio, when com- 
pared to other black hole candidates (e.g. , ICorbel et alj 200< 
20031 iGallo et al.ll200l IFender et aljl2010t ICoriat etal.ll20l1 
The source also appears to be X-ray loud relative to the NIR, 
as in both the hard (epoch 1) and the steep power law (epochs 
2,3) states an extrapolation of the X-ray spectra overestimates 
— significantly in the steep power law state — the NIR flux, 
requiring a low-energy break of the X-ray spectrum. An excess 
X-ray emission is also suggested by a comparison of the NIR 
and X-ray luminosi ties: assuming a distance of 2.6 + 0.7 kpc 
dHoman et all 12006). we find NIR and X-ray luminosities in the 
hard state (epoch 1) of Lnir ~ 10 33 erg s _1 and Lx ~ 10 37 erg s _1 
(each with a 40% error, dominated by the uncertainty on dis- 
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tance). These luminositi es fall under the L MXB (black hole) 
hard state correlation of [Russell e t al. (2006), implying that ei- 
ther the source is X-ray loud (by > 3 orders of magnitude) or 
that the source is NIR dim (by > 2 orders of magnitude). For 
these luminosities to fall on the correlation, the distance would 
have to be a factor of 30, or more, greater (~ 80kpc), which 
is of course unphysical compared to the size of the Galaxy 
(~ 35kpc diameter); indeed, considering the source's position 
in the Galactic plane (/, b = 336.7° -3.4°) any distance much 
greater than 8.5 kpc is unlikely. In the soft state (which has 
luminosities similar to those of the hard state) the source ap- 
pears to be NIR and opti cally dim rela tive to the distribution of 
other soft state LMXBs (Russell et al.l), but not significantly so. 
Furthermore, as noted by iGallo et al.l (120081) . the late-time qui- 
escent X-ray luminosity (Lx ~ 3 x 10 30 ergs _1 ) is at the lower 
end of the distribution of quiescent luminosities for black hole 
binaries, though, when compared with the quiescent optical lu- 
minosity (Lppticai sa 10 30 erg s" 1 ), the X-ray luminosity falls on 
the correlation of iRussell et all albeit at the low end. Obviously, 
because of NIR/optical observational bias (i.e., we only observe 
the brightest sources), caution has to be exercised when implying 
that a source is of a low luminosity compared to a distribution as 
we do here for the soft and quiescent states. However, this obser- 
vational bias should not affect our comparison to the hard state 
correlation. Given that the hard state X-ray luminosity of this 
source is very similar to the range of luminosities presented in 
IRussell et dj (10 35 -10 38 ergs~ 1 ), while the NIR luminosity is at 
the low end of the observed range (at the quiescent level of other 
LMXBs), it seems more likely that the source is NIR, and hence 
radio, dim (a t least quanti t atively ) as opposed to X-ray loud, as 
suggested bv lCorbel et all (120041) . 

While both the NIR and radio fall below their respective 
hard state correlations with the X-ray (in fact, this is the first 
example of an outlier to the NIR/optical correlation), the calcu- 
lated NIR luminosities are more consistent with the correlation 
for the hard states of neutron star LMXBs (Russell e t"aTl l2006h 
— in contrast to the canonical black hole behaviour of the ob- 
served X-ray properties of this source (section[T]i. Likewise, the 
radio luminosities of the source fall on the efficient branch of 
the radio-X-ray correlation and are similar to the luminosities 
of neutron stars (ICoriat et all 1201 lh . The accretion flow in the 
hard state of black hole binaries is generally radiatively ineffi- 
cient and as such a significant amount of the gravitational poten- 
tial may be co nverted to relat ivistic jet flow (i.e., kinematically 
efficient; e.g., IGallo et al.ll2005h . although that energy may also 
pass unobserved across t he black hole's event h orizon, advected 
onto the compact object (Narava n et all fl996T) . Conversely, the 
accretion flow in neutron star binaries (and the soft state of black 
holes) is radiatively efficient and thus thermal accretion emission 
is a more effective method of en ergy radiation than a jet (i.e., 
kinematically inefficient, see e.g., Fe nder et al.ll2.Q03b . Hence, it 
is likely that the detected jet of XTE J 1650 -500 is weak due 
to radiatively efficient accretion (this is what ICoriat et"ai1 refer 
to, qualita tively, as X-ray lou d). This is al so consistent with 
the radio dCorbel et all 12004 and X-ray dRossi et all l2004t 
iMontanari et all 120091; ICorbel etafl 120041) light curves, which 
peak (except at very soft energies, where the soft state thermal 
component is dominant) by the time of our first observation, and 
it is consistent with the X-ray definitions of state since our first 
observation is at the very end of the initial hard state. Thus what 
we observe in both radio and NIR is a weak jet that may be in 
the process of being "turned off" or quenched due to the state 
transition from hard to intermediate as the accretion flow tran- 
sitions from geometrically thick, optically thin and radiatively 



inefficient (but kinematically efficient) to geometrically thin, op- 
tically thick and radiatively efficient. 

5. Conclusions 

The data of XTE J1650-500 suggest that the jet contributes to 
the NIR flux, at least at early times, before being quenched com- 
pletely by the time of the later, steep power law state observation. 
However, we find that the NIR flux falls off much more gradu- 
ally than the radio and maintains a high level of emission even 
as the jet/radio emission goes to zero — suggesting that there is 
also significant contribution from X-ray reprocessing and pos- 
sibly the accretion disk and/or corona, though we cannot distin- 
guish which of these might dominate due to a poorly constrained 
SED and a lack of indicative features of either, such as spectral 
lines or photometric variability. 

The observed radio and NIR luminosities during the hard 
state are orders of magnitude lower than expected for a source 
of this X-ray luminosity, and while several outliers to the radio- 
X-ray correlation of the hard state of LMXBs exist, this is the 
first such example of an outlier to the NIR/optical-X-ray cor- 
relation. The NIR/radio dim emission suggests that the jet of 
XTE J 1650-500 is much weaker than in the hard states of other 
black hole LMXBs, possibly because we observe the jet as it is 
being "turned off", or quenched, at the state transition. 
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